Entanglement of two atomic samples by quantum non-demolition measurements. 
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This paper presents simulations of the state vector dynamics for a pair of atomic samples which 
are being probed by phase shift measurements on an optical beam passing through both samples. 
We show how measurements, which are sensitive to different atomic components, serve to prepare 
states which are close to being maximally entangled. 
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I. INTRODUCTION 

The population of an atomic state |1) can be probed 
by a phase shift measurement on a field that couples 
|1) non-resonantly to another atomic state. This im- 
plies that for an atomic sample with all atoms popu- 
lating two long lived states |1) and |2), it is possible 
to count non-dcstructively the number ni of atoms in 
The state vector or density matrix of the sample 
with, e.g., an initial binomial distribution of populations 
ni will be modified, as the quantum mechanical uncer- 
tainty of the number ni is reduced during the measure- 
ment HJ^L This effect has been demonstrated experi- 
mentally 3I I and further theoretical proposals have 
addressed the possibility to entangle pairs of samples by 
performing joint phase shift measurements, where light 
is propagated through both samples and thus provides 
information about total occupancies of various states. In 
particular it has been possible to prepare EPR-correlated 
samples where the entanglement could be experimen- 
tally proven by a quantitative test involving spin noise 
measurements ||, |9[. 

The theoretical analysis in |^ addresses the behav- 
ior of collective atomic spin operators due to their in- 
teraction with the field operators (in the Heisenberg pic- 
ture). Following the analysis in Q, in the present paper 
we present wave function simulations, in which we dis- 
play how the sequence of photo-detection events gradu- 
ally modifies the state vector of the samples. The two 
approaches are of course equivalent, but they may pro- 
vide different insights, and in addition the present state 
vector approach does not rely on the atomic operators 
being approximated by harmonic oscillators, i.e., we can 
treat cases where the mean polarization changes signif- 
icantly in the samples, which is not practically possible 
in the operator formulation. 

The paper is organized as follows. In Sec. II, we intro- 
duce our physical model for the measurement of optical 



phase shifts, and we present the formal description of the 
state vector evolution conditioned on the photo-detection 
record. In Sec. Ill , we show results of simulations where 



the atoms are first all prepared in superposition states of 
the two internal states prior to detection of the number 
of atoms in state We show that the resulting state 
vector is entangled, and we show that a subsequent spin 
rotation of all atoms followed by new measurements will 
lead to stronger entanglement. In Sec. IV, we show sim- 
ulations where the atomic samples are subject to contin- 
uous spin rotation during measurement, corresponding 
to the experimental situation of Ref. 0. In Sec. ^ we 
present an analysis and interpretation of the results, and 
Sec. VI concludes the paper. 



II. OPTICAL PHASE SHIFT MEASUREMENTS 

In our model we are considering clouds of atoms whose 
relevant level structure consists of three levels: two de- 
generate ground states, |1) and |2), and an exited state, 
|3). We suppose that the only permitted transition is 
|1) |3), excited by the laser beam passing through the 
cloud and followed by spontaneous decay with emission 
rate 7, cf.. Fig. |l|. 
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FIG. 1: Level structure of the atoms. The states |1) and |2) 
are stable states, |1) is coupled off-resonantly by the probe to 
the excited state |3). 
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Our scheme is based on the detection of photons pass- 
ing through the clouds. The scheme only requires classi- 
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cal input fields, since the state reduction following every 
single photodetection event, a posteriori^ extracts the ef- 
fect of the interaction of a single photon with the atomic 
samples. There are different possibilities to create a de- 
tection scheme that measures the population of a single 
atomic state, or the population difference between two 
atomic states. One can use light with two polarization 
components that interact differently with the two internal 
states, and hence the different phase shifts lead to a po- 
larization rotation signal; or, one can use frequency mod- 
ulated light, where one frequency component is closer to 
resonance than the other, and hence a phase difference 
between the two results from the interaction between the 
atoms. For ease of presentation, we follow the descrip- 
tion in [P, and imagine an interferomctric setup, where a 
field is split in two paths, one that interacts with the 
atomic samples and one that propagates through free 
space. This experimental set-up is schematically shown 
in Fig. ||. A phase difference between these two fields due 
to the atoms can be resolved by the intensities measured 
by photo-detectors and D_ in the two output ports of 
the interferometer. We may enclose the atomic samples 
in optical cavities to enhance the interaction with the 
atoms by passing the light through each sample many 
times. This will also reduce the effect of spontaneous 
scattering which will hence be omitted from the present 
analysis. As shown in |^ wave function simulations can 
incorporate spontaneous scattering, and it can also be 
estimated in the operator formulation p]. 



where g is the atom-photon interaction strength, and ol 
is the annihilation operator for the field component in- 
teracting with the atoms. 

If A is the detuning between the light field and the 
atomic transition frequency (E^—Ei)/ h, the second order 
transition amplitude for the interaction process is: 



y 2 



(2) 



where the indices i and / represent the initial and final 
state of the process, which in our case is the same state 
of an atom in 1 1) and one photon of the radiation field. 

In the case of off resonant scattering, with A ^ 7, we 
get 



A' 



(3) 



which is just the energy shift of state |1) due to the ab- 
sorbed and re-emitted photon. Hence, the time evolu- 
tion of the atom-photon state assumes the form |1) (g) 

\nL = 1) -^e-'i^\l)(g)\nL = 1). 
2 

Defining x = ' ''^'^ write an effective Hamiltonian 
describing the interaction between a single mode light 
field and a gas of N atoms: 

Heff ^ J2^\l)n{MalaL ^Y.^a[\a[aL. (4) 





By introducing the atomic spin operators for the 
ground states 



/ n n \ 

Jnz — 2 ^""2.2 ^ "11} 



Jn+ = -^cr^i 



the effective Hamiltonian (W) becomes 



H, 



1 



N 
n—1 
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FIG. 2: Atoms occupying the internal state |1) in the two 
samples interact with the light field which is incident from 
the left in the figure. The phase shift of the light field due 
to interaction with these atoms is registered by the difi'erent 
photo-currents in the two detectors. 

Let us introduce the atomic transition operators afj = 
K)„0|j where i{j) = 1, 2, 3 and n enumerates the atoms. 
In the dipole and rotating wave approximation, the in- 
teraction Hamiltonian is: 



t 



(1) 



In an ensemble of atoms, where each atom is ini- 
tially prepared in the same state, and where the inter- 
action with the surroundings is identically the same for 
all atoms, the collective atomic state retains the full per- 
mutation symmetry, and it is convenient to expand this 
collective state on eigenstates for the effective collective 
angular momentum: 



(9) 



M=-J 



where J 



M 



Y is the total angular momentum, and 



i {n2 — ni) is the eigenvalue of the operator 
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X^n jriz- Collective raising and lowering operators, and 
the corresponding cartesian x- and y- components of the 
collective angular momentum, are defined as similar sums 
over all atoms in the sample. 

A photon incident on the interferometer splits as a su- 



perposition of a photon in the upper path and a photon 
in the lower path, -^{\U)pf^ + \L)pf^), and the time evolu- 
tion of the collective atomic state due to the interaction 
with this photon for a duration r is given by 



Mr)) 



ph-\-at 



M=-,J 



I 



(10) 



where (f - M) = (iii±ii2 „ "2^) ^ Hence, as 
claimed above, the phase shift of the single photon state 
\L)ph is proportional to the total population of state 



A. Entanglement due to measurement 

In 1^ the interaction with a single sample is analyzed in 
detail, and the possibility to achieve spin squeezing due 
to the interaction with the atoms is analyzed. We now 
turn our attention to a system of two atomic ensembles 
of Ni and N2 atoms, respectively. Within each sample 
we assume permutation symmetry among the atoms, i.e., 
each sample is represented by the collective spin states. 



J 



Mi,M2 



\u) 



introduced above. 

The state vector of the samples can be expanded on 
product state wave functions of the two ensembles 

|*>= J2 ^M„M2\Ml,M2), (11) 
Mi,Af2 

where \Mi,M2) = \ Ji = iVi/2,Mi) ® IJ2 = A^2/2,M2), 
and where we assume an initial product state Ami,M2 = 
AmiAm2 before we submit the system to the incident 
field. The photon interacts with each of the ensembles 
in the lower arm of the interferometer, and the state of 
the system after the interaction, but before the photo- 
detection, is 



ph 



-iMi+M2)]xT 



ph 



V2 



\Mi,M2 



(12) 



The detectors D+ and I?_ click, when they detect a pho- a+ and a_. It is thus useful to rewrite the photon-atom 
ton in one of the states |±)„,, = -j= (\U).„,^ ± \L)A , and state exphcitly in terms of the two components distin- 

. , , . , , . , , . , 1 , • , guished by the photodetectors: 

associated with this detection event the photon is de- 
stroyed by the action of one of the annihilation operators 

I 



l*(T')>p/i+at = E -^Afi3/2 ( - 
Mi,M2 \ 



Mi,M2 



|Mi,M2)<»|+)p^ 



/ I ^ e~i^^-iM,+M2)]xr\ 
.M2 ^ |Mi,M2) 



/ph ■ 



(13) 



The probabilities that the photon is detected in state \+)pk or in state are 
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i\/l,J\/2 
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(14) 



(15) 



Thus, after the detection of a photon, the state of the 
two ensembles of atoms is projected, with probabihty 7r_|_ 
and 7r_ respectively, into one of the following states 



1 



E 



M2 



1 + e 
1 



Mi,M2 
_ rjvi+jva 



1 -e 



E 



-(Afi+Af2)]xr^ |Mi,Af2) (16) 

-4j\/l,Af2 



Ml, Ah 



iMi+M2)]xr\ |Mi,Af2),(17) 



where C± are normalization factors. From Eqs. ( p^ ) 
and (l^' we observe the entanglement between the two 
atomic ensembles emerge. 

The detection procedure, and the corresponding wave 
function updating, is repeated for a number Nph of pho- 
tons, being detected one at a time. Hence, by defining 
the following entangling factors 



T±{Mi+M2) 



lie" 



•^-{Mi+M2)]xr 



, (18) 



after Nph detected photons, of which Af+ are detected by 
D+ and = Nph — N+ are detected by _D_ , the state 
vector of the samples is 



^ ^M,^M2[^+(Afl+M2)]^+ 



Mi,M2 



X [T^{Mi + M2)]^^'"^~^+'> \Mi,M2) . (19) 

The square norm of the entangling factor in equation 
® is 

r /w_A/r,„ \i2^+ 

B{A'h2) = 



N ~ M12 
cos I XT 



N - M12 
sm I XT 



-,2{Nj,H-N+) 



(20) 



where Afi2 = Mi + Af2, and where, for simplicity, we 
assume the same number of atoms in each sample — 
Ni — N2- For large Nph this function is very peaked with 
maxima at values M12 = M obeying 



,21, 



Equation (|T|) has multiple solutions for M, due to both 
the periodicity of the ian-function and the double sign in 
the right hand side, but we can remove this ambiguity in 
the solutions if we ensure that only a single maximum is 
compatible with the initial distribution of the two sam- 
ples in Ml and -M2. Calculating the second derivative of 
B{Mi2) around the maximum values, we get an estimate 
of the width of the peak in M12 



dl,^^B{Mi2)\ 



Mi2 = M 



-{xrfNphB{Mi2) 



(22) 



If the number of atoms is large enough, we can approxi- 
mate B{Mi2) with a Gaussian, and with this approxima- 
tion we obtain the r.m.s width 



1 



(23) 



Thus, the more photons that are detected, the more 
peaked is the distribution in -M12. 



III. NUMERICAL SIMULATION, 
CONSECUTIVE MEASUREMENTS OF J^i + ,7^2 

AND Jyl - Jy2 

In this section, the detection model described in the 
previous section will be implemented in a numerical sim- 
ulation in which the two samples have the same, not 
very large, number N of atoms and the initial state is 
the eigenstate of the operator Jxi + Jx2 with eigenvalue 
Ji + J2 = 2J = N, i.e., all the atoms are initially pre- 
pared in the superposition state -^(|1) + |2)). The cor- 
responding amplitudes in the basis eigcnstatcs of the z- 
component of the collective angular momentum operators 
|J, Ali) are 



= 

where i — 1,2. 
atoms in state 

the amplitude (|2|) is \A 




(24) 



Expressed in terms of the number of 
|1) in each sample, n\, the square of 



N 



means that in each ensemble the atoms are distributed 
in state |1) and |2) according to a binomial distribution 
with probability i. 

Beginning with an initial state characterized by am- 
plitudes (p4h, we proceed with a series of photo detec- 
tions that project the state of the system according to 
Eqs. ( p^ and (17). After a large number of photons 
Nph has been detected, the uncertainty in Afi2 is re- 
duced, i.e., we are almost in one of the eigenstates of 
Jiz + J2z- Fig- m shows the behavior of the variance 
A^(Ji2 + J2z) as a function of the number of detected 
photons for three different simulation records. It is evi- 
dent that when the number of detected photons is large 
enough, A^( Ji^ + J2z) goes to zero, i.e we have a definite 
value of M = Mi + A^2- In the simulations we assume a 
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phase angle = 0.24. Such a large phase shift on the 
atomic state due to interaction with a single photon is not 
realistic in experiments with freely propagating fields. In 
our simulations a smaller value of xr just implies that 
more photons have to be detected to achieve the same 
reduction in A^(Ji2 + J2z), c.f., (p3|). 
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FIG. 3: Variance of Jiz + J2z as function of the number of 
detected photons. Each sample consists of 20 atoms in the 
simulations, and the phase angle xt per atom and photon is 
0.24 in the simulation. 

To quantify the entanglement between the two sam- 
ples we apply the definition by Bennett et al of the 
entanglement between two systems in a pure state I^E"): 



8 = Tr{pi log2 pi) = Tr{p2 loga P2), 



(25) 



where pi = Tr2 \^) (vE'| is the reduced density matrix of 
system 1. (Similarly for ^2). If is the number of atoms 
in each sample, and we restrict ourselves to states which 
are symmetric under permutations inside samples, the 
quantity £ takes values between zero for a product state, 
and log2(iV + 1) for a maximally entangled state of the 
samples. 

Fig. ^ shows the behavior of the entropy in three dif- 
ferent simulations. The left hand side of each figure, for 
Nph < 500 corresponds to the detection, analyzed above. 
To prove that the samples are entangled and not just clas- 
sically correlated, in the experiments of Q a second ob- 
servable was introduced, for which the quantum mechani- 
cal uncertainty was similarly reduced, and measurements 
on both sets of observables conclusively demonstrate the 
entanglement H, §|. In our theoretical calculation, we 
know of course already that the state is entangled, but 
as we shall see, we can increase the entanglement by tak- 
ing a second round of measurements. 

Thus, after the detection of 500 photons, we apply op- 
posite rotations to the atomic samples, and we proceed 
with similar measurements as before, which effectively 
J2y in the non-rotated frame. After the 



rotations, the full permutation symmetry of the system is 
broken, i.e the total angular momentum = J^+Jy+J^ 
is not conserved. The final state of the two samples is 
still an eigenstate of both and J| , but with different 
values of the total J. Our formalism (25,26) is already 
prepared to handle that situation, since the basis vectors 
in the expansion are kept on the form of product states 
rather than angular momentum coupled states of the two 
samples. 

As illustrated by the right hand parts of the panels in 
Fig. ^, the entanglement had essentially saturated after 
the first 500 detection events, but after the rotation, it 
grows again to saturate, typically, at a higher level. In 
Fig. ^(a) we show an example of the evolution of the 
entropy in which the final value {Enj,^ = 4.3193) is very 
close to that of the maximally entangled state (log2(20-|- 
1) = 4.3923). In Fig. |(c), however, we show a case in 
which the final entropy value is reduced by the second 
set of rotations. In Sec. |^, we shall return to an analysis 
of these findings. Finally, Fig. ^(d) presents the average 
of the entropy over 50 simulations. 
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FIG. 4: Entanglement of two samples of each A'^ = 20 
atoms. After detection of the first A'p^ = 500 photons, the 
samples are rotated ±90 degrees in spin space, and the sub- 
sequent detection events lead, typically, to further entangle- 
ment. Figs. 4 (a-c) show results of three different simulation 
records. Fig. 4(d) shows the average over 50 such simulations. 
XT = 0.24 as in Fig. 3. 



IV. MEASUREMENTS WITH CONTINUOUS 
ROTATION OF SPIN COMPONENTS 

As we just saw that measurements of two different sets 
of operators typically lead to an increase of the amount 
of entanglement, it is natural to consider a continuous 
change between these operators. In fact, in the experi- 
mental work by Julsgaard et al [Q , such a continuous ro- 
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tation was induced for purely practical reasons (so that 
the relevant signal could be picked up in the absence of 
technical noise at a high frequency). In this section we 
present the results of simulations in which continuous op- 
posite rotations are applied to the atomic spins of the two 



samples as the photo detection proceeds. If ^jv/^^f^ 
wave function amplitude of the system after Nph photo 
detections, the updated state vector rotated by the angle 
±9 prior to the subsequent detection is 



E 

m',,m!, 



y v'^ , (9) v'- 



Mi,M2 



(26) 



where D^^^^, (9) = (M„ J,| e-*^"(«) 



Ji,M^j is the ro- 
tation matrix clement for sample i — 1,2. The wave 
function updating algorithm has the same structure as 
described in section II. In a realistic experimental situ- 
ation, one would detect photons according to a Poisson 
process, and with a constant rotation frequency, induced, 
e.g., by application of opposite DC magnetic fields onto 
the atoms, this would lead to small rotation angles with 
an exponential distribution law. For simplicity, however, 
we apply the same small rotation angle 9 between each 
detection event. 



In Fig. H we show realizations of the detection scheme 
with small rotations. The rotation angle between subse- 
quent photodetection events is 9 — ^. Fig. ^(a) shows 
a case where the entropy evolves to the maximum value: 
Sff^^ = 4.3923. In Fig. |(b) and |(c) we present 'typi- 
cal' and 'worst case' results of the simulations, while in 
Fig. ||(d) we provide the average over 50 simulations. 

We observe in these figure that compared to the results 
shown in the previous section, there is a faster increase 
of the entanglement towards an almost constant level, 
which varies from one simulation to the other. 
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FIG. 5; Entanglement of atomic samples with 20 atoms. The 
spins are rotated in opposite directions around the a;-axis by 
the angle = j after each detection event. Insets (a) and 
(b) show the evolution of the overlap | ('I'ol ^')iVphP between 
the state of the samples and the maximally entangled state 
(See section (0)). In (a) 1*)^^ ^ converges towards |\['o) and 

further detections have no effect on the state vector. In (b) 
the state gradually gets rid of its component along |^o), and 
it subsequently evolves in the orthogonal subspace of |^o). 
In (c) is shown a record where the entanglement is smaller 
than in the simulations in (a) and (b), and in (d) is shown 
the average over 50 simulations. 



V. ANALYSIS AND INTERPRETATION 

In order to interpret the results obtained in the pre- 
vious sections, let us first address the achievements of 
probing two operators rather than a single one. The im- 
provement is quite easily understood, if one notes that 
the detection of Jzi + Jz2 effectively produces an approx- 
imate eigenstate of this operator. This state has am- 
plitudes on different Mi and M2 states with Mi + M2 
fixed by the measured value, but the eigenspace is de- 
generate, and the amplitudes are simply proportional to 
the ones in the initial state. Due to the initial bino- 
mial distribution on Mi and M2, the distribution over, 
e.g., Ml will therefore have a width of approximately 
ViV. The reduced density matrix has the correspond- 
ing number of non-vanishing populations, suggestive of 
£ ^ log2 -y/ZV = 0.51og2iV, half of the maximal value. 
This argument accounts for the first plateau reached in 
Fig. 4. 

The measurement of the z-components causes a broad- 
ening of the distribution on Jy eigenstates beyond the ini- 
tial distribution, which was also binomial in that basis. 
The subsequent measurement of Jyi — Jy2 will produce 
a state with Miy — M2y fixed by the measurement, but 
within the degenerate space of states with this fixed value 
the distribution on Miy of the reduced density matrix is 
broader than -s/iV, and the entanglement is correspond- 
ingly larger. This accounts for the increase of entangle- 
ment seen in most simulations. 

For large collections of atoms with large mean values 
of the collective operators, the orthogonal spin com- 
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ponents Jz and Jy are well approximated by effective 
position and momentum operators, and for two particles, 
the pair of combinations of position and momentum 
operators xi + X2 and pi — p2 commute, i.e., they can 
both be measured with high precision. But, this is 
not an exact replacement, and in our simulations with 
fewer atoms, we see significant deviations from this 
picture. The commutator of Jzi + Jz2 and Jyi — Jy2 is 
proportional to the operator J^i — Jx2 which does not 
vanish, and in general, measurements sensitive to one 
of these operators are complimentary to measurements 
sensitive to the other one. It is possible, however, to find 
a single ]omi eigenstate of the operators 



(Jlx- J2r.)|*0) = 0, 

(Jly- J2,)|*0) = 0. 

Rewriting these equations as 

(Jl+- J2+)|vI'o) = 0, 

(Jl- - J2-)|*0> - 0, 



(27) 
(28) 



(29) 
(30) 



it is easy to check that the solution is the maximally 
correlated state 



^/2J 



1 

= V |M,-M) 



(31) 



M=-J 



Note that |^'o) also satisfies (Jiz + J2z)|^'o) = 0, and in 
fact, all spin components have vanishing mean values in 
this state. 

When we simulate the detection of phase shifts propor- 
tional to J1Z + J2Z and Jiy — J2y , or combinations of these 
operators, there is a chance, that the state vector is grad- 
ually projected onto |5'o), and this state is unaffected by 
all future photodetection events. Indeed, in some of our 
simulation records, we precisely see the robust generation 
of the maximally entangled state, cf.. Fig. ||(a). 

If one has not collapsed into j^E'o) after a large num- 
ber of photons have been detected, the state vector is 
instead orthogonal to that state, and one will never ar- 
rive at the maximally entangled state. The insets of the 
figures show the wave function overlap with |5'o), and 
in all simulations this quantity converges to either unity, 
as in Fig. ||(a), or to zero, as in Fig. ^(b). The overlap 
between 1^*0) and our initial state suggests that in one 
out of iV -|- 1 realizations of the experiment, one should 
produce that particular state, and this is confirmed by 
our simulations. 

There are no other joint eigenstates of the pair of col- 
lective operators, and hence state vectors orthogonal to 
l^'o) do not evolve into any specific state: measurements 
on the system have different outputs which affect the 
state vector in different ways. We do note, however, that 
fairly strong entanglement is observed in many realiza- 
tions, and that this entanglement is almost constant over 
time. This is suggestive of families of states with rela- 
tively well defined values of the operators Jiz + J2Z and 



Jiy — J2y, which is allowed by Heisenberg's uncertainty 
relation as long as the expectation value of the operator 
Jix ~ J2X is small. Figs. 6-8 present the three mean val- 
ues of Jix — J2X, Jiy — J2y and Jiz -\- J2Z for the evolution 
leading to the maximally entangled state, a very entan- 
gled state and a less entangled state studied in Fig. ^ 
(a-c). The picture is precisely as expected: the maxi- 
mally entangled state does not change with time and the 
mean values vanish forever, a strongly entangled non- 
stationary state performs almost regular oscillations in a 
limited part of Hilbert space restricting the mean values 
to a similar small oscillatory behaviour, and states with 
little entanglement show more dramatic time dependence 
and the mean values explore a wide range of values. It 
was unexpected that the non-maximally entangled states 
seem to be consistently much entangled or little entan- 
gled for long detection sequences. Further studies of this 
dynamics would be very interesting. 
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FIG. 6: Evolution of (Ji^ — J22,), (Jiy — J2y) and (Jiz -f- J2z) 
corresponding to the simulation shown in Fig. ^(a). 



VI. CONCLUSION 

We have presented a quantitative wave function analy- 
sis of the entanglement created by total population mea- 
surements on separate atomic samples by means of op- 
tical phase shifts. The analysis included wave function 
simulations where small samples were exposed to inter- 
action with the optical fields in an interferometric set-up. 
We recall that more realistic experimental set-ups can be 
made , but they will be treated by the same formalism 
and lead to the same results. Our choice of parameters 
(only 20 atoms, high atom-field coupling) was not meant 
to correspond to a specific experiment, but it brings out 
results, that we may translate also to realistic experimen- 
tal regimes. Note, however, that the present work should 
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FIG. 7: Evolution of (Ji^ — J2x), {J\y — Jiy) and (Ji^ + J^z) 
corresponding to the simulation shown in Fig. |(b). 

be extended to include also spontaneous emission by the 
atoms, if one wants to model experiments, where this has 
significant effects on the dynamics. 

The entanglement protocol has already been imple- 
mented experimentally , and a theoretical analysis un- 
der an harmonic oscillator approximation has been pre- 
sented. On one side our analysis supplements this exist- 
ing analysis with a state vector perspective, and on the 
other hand it provides an analysis valid beyond the oscil- 
lator approximation, where in particular the emergence 



of a single maximally entangled state, and families of 
non-stationary states with different degrees of entangle- 
ment were identified. Since the experimentalist in princi- 
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FIG. 8: Evolution of (Jix — J^x), {J\y — J2y) and (Ji^ + J2z) 
corresponding to the simulation shown in Fig. ^(c). 

pie knows the state vector conditioned on the outcome of 
the detection, it seems interesting to introduce feed-back 
I p^ , either by just resetting the system and start over 
again if an only weakly entangled state is prepared, or 
by suitably hitting the system during measurements. 
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